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Abstract 

The differential cross sections for low-energy muonic hydrogen atom scattering in solid molecular 
H2, D2 and T2 targets under low pressure have been calculated for various temperatures. The 
poly crystalline fee and hep structure of the solid hydrogenic targets are considered. The Bragg 
and phonon scattering processes are described using the Debye model of a solid. The calculated 
cross sections are used for Monte Carlo simulations of the muonic atom slowing down in these 
targets. They have been successfully applied for a description of the production of the muonic 
atom beams in the multilayer hydrogenic crystals. 
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I. INTRODUCTION 



The aim of this paper is to calculate the differential cross for muonic hydrogen atom scat- 
tering in low-pressure molecular hydrogenic crystals and to study the muonic atom deceler- 
ation in such crystals. The solid hydrogenic targets have been used for the production of the 
muonic hydrogen atom beams, which have been applied to the time-of-flight measurements 
of various muonic-atom and muonic-molecular processes (see e.g., Refs. H, H, 0, H, 0]). 
The bulk hydrogenic crystals are employed in the investigations of muon- catalyzed fusion of 
the hydrogen isotopes at high target densities and for various populations of the rotational 
molecular levels [J @, 0, [lO, 11, 111, 13 1. Hydrogenic solid targets are also used in a novel 



method of spectroscopy of the radioactive muonic atoms 14. isl . 

Monte Carlo simulations of such experiments, which used the perfect-gas model of a tar- 
get and the available differential cross sections 16, l3| for the muonic atom scattering 
from hydrogen-isotope nuclei, did not reproduce the experimental low-energy (< 1 eV) 
data. Moreover, even the use of the cross sections for scattering from isolated hydrogenic 
molecules [isl, [11] did not lead to an agreement between theory and experiment Thus, 
it is necessary to take into account condensed-state effects of the muonic atom scattering in 
the solids. In particular, interpreting the experimental data and planning new experiments 
demand the knowledge of the differential cross sections for scattering in the hydrogenic 
crystals. 

A few examples of the total cross sections for pfi and dn atom scattering in solid H2 and 
D2 at a fixed target temperature T = 3 K were shown in Ref. [23]. Below, a method of 
calculation of the partial differential cross sections for muonic hydrogen atom scattering in 
the polycrystalline hydrogenic targets at various temperatures is presented in detail. These 
cross sections are evaluated using the amplitudes for muonic atom scattering from isolated 
hydrogenic molecules [l9| and the Van Hove response function S [2l[. The isotropic Debye 
model of the solid is employed. 

In Sec. [ni the coherent and incoherent cross sections for the scattering in hydrogenic 
crystals are expressed by the scattering amplitudes calculated for the isolated molecules. The 
incoherent cross section, which includes the elastic scattering, incoherent phonon scattering, 
and rotational-vibrational transitions in the target molecules is discussed in Sec. IIIII The 
Bragg scattering for the fee and hep structures, which are observed in hydrogenic crystals, 
is considered in Sec. IIVI The inelastic coherent cross section, which leads to the creation 
or annihilation of one phonon in a muonic atom collision with the crystal is calculated 
in Sec. |Vl Some examples of the calculated cross sections for the homogeneous hydrogenic 
targets at various temperatures are shown in Sec. IVI[ Most examples are given for the so- 
called "normal" targets nH2, nD2, and nT2, which are characterized by the statistical 
distribution of the molecular rotational levels K = and 1. 

Symmetry of the wave function of a homonuclear hydrogenic molecule is definite. In 
the hydrogen case, the total nuclear spin I of H2 equals or 1. The singlet spin state 
/ = is asymmetric. As a result, only the symmetric spatial wave functions with even 
values of the rotational number K are allowed. The H2 molecules in such states are called 
parahydrogen molecules. For the symmetric triplet state 1 = 1, only odd values of K 
can occur. This defines the orthohydrogen molecules. The situation is similar for the T2 
molecules, since the spin s of proton or triton is equal to 1/2. In the deuterium case, one 
has s = 1. The total two-deuteron wave function must be symmetric. For the symmetric 
spin states / = 0, 2, the spatial wave function is symmetric, which corresponds to even K. 



2 



The D2 molecules in such states are called orthodeuterium molecules. The paradeuterium 
molecules are characterized by / = 1 and odd K. The ortho and parastates of the hydrogenic 
molecules are remarkably stable, in the absence of a catalyst ^]. As a result, during rapid 
cooling of the equilibrated nH2, nD2, or nT2 gas, the even-i^ states deexcite to K = 
and the odd-i^ states deexcite to = 1. The ortho-para transitions can occur in certain 
collisions with muonic hydrogen atoms. This is possible if the muonic atom interaction with 
a homonuclear molecule is spin-dependent. Such a situation takes place when the hydrogen 
isotope in the impinging atom is identical with those in the homonuclear target molecule, 
e.g., in the p/x + H2 scattering [l6l, 13]. In this case, the exchange forces between the three 



identical nuclei should be taken into account. In particular, they lead to a high probability 
of the muon exchange between two nuclei taking part in direct collision. When the spin 
projections of these nuclei are opposite, the nuclear spin / of the target molecule can change 
in the collision process. As a result, the simultaneous rotational ortho-para transition takes 
place. Such spin effects are included in the presented calculations. 

The obtained differential cross sections are used for Monte Carlo simulations of the muonic 
atom deceleration. Some results of these simulations are presented in Sec. I VI II 



II. COHERENT AND INCOHERENT CROSS SECTIONS 

A muonic hydrogen atom a/i can be approximately treated as a small neutron-like parti- 
cle. Therefore, the methods derived for the description of neutron scattering in condensed 
matter can be applied to a/i scattering in dense hydrogenic targets. Below, the Van Hove 
formalism [2l| is adapted to the calculation of the differential cross sections for a/i scattering 
in the hydrogenic crystals. In this formalism, the cross sections are expressed in terms of the 
response function iS, which depends solely on properties of a given target for fixed momen- 
tum and energy transfers. A definition of this function involves both quantum-mechanical 
and statistical averaging over the target states at temperature T. 

It is assumed that there is no coupling between the translational and collective motions 
of the molecules in the condensed target and the molecular rotations and vibrations. The 
internal molecular degrees of freedom are already included in a single-molecule process 1^ 



and, therefore, do not enter into the response function. It is assumed that a bulk dense 
target is kept at a sufficiently low pressure, so that distortions of a single bound molecule 
due to the interactions with neighbors can be neglected. This assumption is fulfilled in 



the case of low-pressure 10 kbar) solid hydrogenic targets [22|, |23|. The mean distance 
between the neighboring molecules in such targets is several times greater than the diameter 
of these molecules. The Van der Waals force between the molecules is weak. As a result, 
the rotational and vibrational numbers remain good quantum numbers, although small 



broadening of certain excited molecular levels takes place [2J] . This broadening is not taken 
into account in the presented calculations. 

Low-pressure hydrogenic solids are quantum molecular crystals, which are characterized 
by a large amplitude of the zero-point vibrations of the molecules in the lattice. The stan- 
dard lattice dynamics can be applied to these crystals, after certain renormalization of the 



molecule-interaction potential |22l . l23l |. Also, the Debye model of a solid can be used as 
a reasonable approximation. 

The wavelength of a very slow (< 10 meV ) rn uonic hydrogen atom is comparable to the 



nearest-molecule distance of about 3.5 A |22|, |23|. Therefore, strong interference effects can 



be observed at such energies. These effects are described using a conventional separation of 
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the total differential cross sections on the incoherent and coherent fractions. The coherent 
scattering takes place only if specific geometrical conditions are fulfilled. Analogous to 
the neutron coherent scattering the coherent cross section for a/x scattering in a solid 
single-isotope hydrogenic target can be written down in the form: 

/ \ _ k' 



dnde 



\ k' 

-\ =Nraol-rO-cohS{K,Uj), (l) 
/ coh 



where A^moi is the number of molecules in the target. The energy transfer u and the mo- 
mentum transfer k to the lattice are, respectively, equal to 



UJ 



e-e' -AE, k = k - k', (2) 



where e and e' denote the initial and final kinetic energies of the scattered muonic atom 
and AE is the sum of the internal-energy changes of a/i and of the target molecules. Vectors k 
and k' stand for the initial and final momenta of a/x. These momenta and collision energies 
are connected by the relations 

e = '^e/Ma, , e' = \k'^/M,^ , (3) 

in which M^^ denotes the mass. The function a^oh in Eq- dl]) is expressed by the ampli- 
tude jF™°^ for a/i scattering from an isolated molecule [l9| 



^coh = . (4) 

The horizontal bar stands here for averaging over a random distribution of the total spin J 
of the a/i+molecule system and over a distribution of the initial rotational states of the 
molecules. It is assumed that there is no correlation between the direction of the molecular 
spin I and the lattice site. 

Incoherent scattering does not include interference effects from different molecules in the 



lattice. The incoherent cross section takes a general form [25 



( a^a,/ l = ^mol T: C^inc Si{hL, u) , (5) 



dQde'),^^ -k 
where 

C^inc = |J^"^°1|^ - l^^l^ (6) 

and the incoherent response function Si{K,uj) is a fraction of the total response func- 
tion S{k,,u!). In the limit of large momentum transfers, the coherent processes disappear, 
so that S{k,,u!) ^ Si{K,uj). The total differential cross section d'^a / dQde' is a sum of the 
coherent ([1]) and incoherent ([5]) cross sections. 

The inelastic scattering processes which change the internal state of a/x or that of the 
target molecule (such as the spin-flip and isotopic-exchange reactions and the rotational- 
vibrational transitions) are fully incoherent processes. No averaging over the states of the 
different target molecules is performed and, therefore, in this case (Jinc reduces to the single- 
molecule squared amplitude 

o"mc = and (Jcoh = . 



4 



When the states of the muonic atom and molecule are not changed during collision, 
(Tcoh and (Jinc are equal to the coherent and incoherent fractions of the elastic cross section 
for a/i scattering from a single molecule. Their values depend on a given choice of the 
hydrogen isotopes, the total spin F of afi, the population of the molecular rotational levels, 
and the collision energy. In particular, when JF™°' does not depend on the spin J' and only 
one rotational state is populated, the scattering is fully coherent 



O'coh 



mol 



■moll 



and 



0, 



For example, such a situation takes place in the case of elastic scattering dfi+B.2 in the 
ground rotational state 7^ = of the H2 molecule. 

In general, both (Xcoh and o"inc can have nonzero values. In Table [H these functions are 
shown in the limit e — > 0, for the cases pfi + B.2, + T2, and dfi + D2 and for several values 
of F and K. 



TABLE I: Coherent 


(Tcoh and incoherent 


(Tine elastic 


cross sections (in 10 '^^ cm^/sr) 


for a single 


hydrogenic molecule 


at the collision ener 


gy e ^ 0. 






process 


F 


K 


0"coh 


•^inc 










88.5 


0.0 







1 


88.7 


0.0 




1 





126.2 


0.0 




1 


1 


124.8 


87.6 










325.7 


0.0 


t^+Ts 





1 


324.9 


0.0 




1 





342.9 


0.0 




1 


1 


342.4 


10.5 




1/2 





64.7 


0.15 




1/2 


1 


64.4 


0.06 




3/2 





64.5 


0.38 




3/2 


1 


64.3 


0.15 



III. INCOHERENT SCATTERING 



The incoherent response function Si can be rigorously calculated for harmonic crystals [21|, 
26( 1 . Using the so-called phonon expansion of Si for a cubic Bravais lattice (one molecule per 



lattice cell), the incoherent cross section ^ takes the form: 



X 



{2wy 



(7) 



n=l 
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where exp(— 21^) denotes the Debye- Waller factor, which is familiar in the theory of neutron 
scattering. This expansion is also a fair approximation for any cubic lattice, e.g., the fee 
structure. The exponent of the Debye- Waller factor is by definition equal to 

2W{k)^{{k.u)')^, (8) 

in which u is the displacement of the molecule from its lattice site and (. . .)t denotes the 
quantum-mechanical and statistical averaging at temperature T. In particular, for a cubic- 
crystal structure 

2W=1{u^)tk\ (9) 

It follows from that 2W does not disappear at T = because the mean square dis- 
placement tends in this limit to a finite value determined by the zero-point vibrations of the 
molecule in the lattice. The functions gn in Eq. ([7]) are defined as 

/X 1 Z(w) . , , 



gn{w) = dw' gi{w - w') gn-i{w') , (10) 

J — oo 

dwgn{w) = 1 , 



where M^^i is the mass of the molecule. The normalized density of vibrational states Z{w) 
in the isotropic Debye model has the following form: 

10 ii w > wj) , 

where wd = k^Oi) is the Debye energy corresponding to the Debye temperature Oj) 
(/cb stands for the Boltzmann constant). For the low-pressure hydrogenic crystals, 0-d is 
on the order of 100 K. Although differences between the masses of the hydrogen isotopes 
are quite large, isotopic effects in for these crystals are quite small, due to quantum 



effects [22|, 1231]. The function nB{w) in Eqs. ( |TOl) denotes the Bose factor 

nB{w) = [exp{w/kBT)-l]-\ (12) 

which determines the phonon population at a given T. 

Note that there is no momentum-conservation condition in the incoherent cross sec- 
tion ([7]). The 6 function in Eq. ([7]) represents the recoil- less scattering from the rigid lattice. 
In the case of elastic scattering {AE = 0), a/i cannot change its energy in this process. 
This is connected with a large mass of the solid target. The next terms give broad dis- 
tributions corresponding to the subsequent multiphonon processes. In particular, the term 
with n = 1 describes incoherent afi scattering with simultaneous creation or annihilation 
of one phonon. The amplitudes of all the processes are proportional to the Debye- Waller 
factor, which decreases with the rising momentum transfer. 

Elastic a/i scattering from a free molecule is described by a single function. In a solid, 
this function is replaced by a set of functions corresponding to the strictly elastic and 
multiphonon processes which are proportional to k^". The same conclusion can be drawn 
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for any incoherent process, e.g., for a rotational excitation of the target molecule. In such 
a scattering event, the nonphonon cross section is connected with the energy transfer AE 
to the molecule. The higher energy transfers are due to the rotational excitation with 
simultaneous phonon creation. 

When the momentum transfer is small {2W ^ 1), only the several lowest terms in Eq. ([7]) 
are important. On the other hand, in the limit 2W ^ 1 (weak binding), many multiphonon 
terms give comparable contributions to the cross section. For sufficiently large k^, it is 
convenient to use the impulse approximation in which Si takes a general Gaussian form 
for any target described by a time- independent Hamiltonian In this approach, the 

incoherent cross section is equal to 



(■ 



\df2de'J. k 

(13) 



inc 



1 

X := exp 

Arv^ 



, 2 



A 



R 



where 



Ar = v/8^t^r/3 (14) 

is the Doppler width of the asymptotic form of Si. The mean kinetic energy of a single 
molecule in the lattice at temperature T is denoted by and is the recoil energy 



For a harmonic solid 



wr = i/c'/M^oi . (15) 

/•OO 

= I / dw Z{w) w [riBiw] + i] . (16) 
Jo 

This energy contains a contribution from the zero-point vibrations of the lattice molecules. 
In the case of low-pressure hydrogenic crystal, one has T /0d <^ 1. In this limit, St and the 
effective target temperature T^s 

Teff=KTAB. (17) 

are well approximated by 

c^T = ^wo, and Teff=f0D, (18) 



respectively. In particular, for a 3-K zero-pressure solid deuterium target with = 108 K, 
one obtains St = 5.2 meV and Tes ~ 40 K. Thus, T^s ^ T, which means that the cross 
sections calculated for 3-K perfect-gas D2 should not be used for a description of the solid D2 
target at the same temperature. Relatively high values of St in the solid H2 and D? targets 



have been experimentally confirmed using deep inelastic neutron scattering [281. l29l. |30|. 

At very high collision energies e ^ wd and e ^ AE, the approximation e' ^ e [k' k) is 
valid. As a result, the cross section f|T3l) averaged over e' tends to the static approximation 



(19) 



This cross section no longer depends on the target structure and is equal to the sum of the 
free-molecule incoherent cross sections. Since, at high energies, coherent effects disappear, 
the total differential cross section dcr/di? is given by Eq. (|T9l) with (Tine replaced by |jF™°^p. 
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The influence of tlie lattice binding on the incoherent scattering is described by the 
following function: 



mol 



M, 



M 



dQde 



(20) 



where M. denotes the reduced mass of the a/x+molecule system. It is assumed here that 
jrraoi jg constant. The function Cine for elastic {^E = 0) pfi scattering in 3-K solid H2 at zero 
pressure is presented in Fig. [T] The different curves show contributions from the nonphonon, 
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FIG. 1: (Color online) Function Cine versus energy e for incoherent elastic scattering of p^u-atom 
in 3-K solid- hydrogen. Subsequent phonon-creation processes are denoted by (one-phonon), 
"+2" (two-phonon) , and "+3" (three-phonon). Label "—1" stands for annihilation of one phonon. 



phonon creation, and phonon annihilation processes. Every single process falls exponentially 
for sufficiently large energies, which is due to the Debye- Waller factor in Eq. ([7]). At energies 
below a few meV, the nonphonon elastic scattering is dominant. The phonon annihilation, 
which is strongly suppressed by the Bose factor at 3 K, prevails over the phonon creation only 
at e < 0.7 meV. The subsequent phonon-creation processes appear when e rises. Finally, 
they dominate the scattering at e ^ wb, where the static-approximation limit is reached 
and Cine = 1- For e ^ , one has Cine = {Ma^/MY = 2.24. In Fig. [21 contributions of the 
nonphonon and phonon fractions of Si to the total cross section are plotted for the rotational 
deexcitation = 1 — of a D2 molecule bound in 3-K solid deuterium, in collision with 
a dfx^F = 3/2) atom. The cross section is proportional to at £ — >■ 0. Even at the 

lowest energies, the phonon processes are more important than in the elastic case (cf. Fig. [T]) 
because the rotational energy release of 7.5 meV is comparable to wd ~ 9 meV. 



IV. COHERENT ELASTIC SCATTERING 

The fraction of the response function S that describes coherent elastic scattering (Bragg 
scattering) is well known in the neutron theory (see e.g., Ref. |,25!]). When applied to Eq. ([T]), 
it leads to the following coherent elastic cross section for a perfect-crystal lattice with Nd 
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FIG. 2: (Color online) Phonon and nonphonon contributions to the cross section for the rotational 
deexcitation = 3/2)+D2(i<" = 1) ^ dfi{F = 3/2) +D2(i^ = 0) in solid deuterium at T = 3 K. 



identical molecules per unit cell of volume Vq: 



X 6{k, — r) exp [— 2iy(K 
Fn{t) being the unit-cell structure factor 



I ' ° 2 

V (r) I ^ = (Tcoh 5Z '^^'^ ■ ^) • 



d=l 



The summation in Eq. fl2Tl) is performed over the reciprocal-lattice vectors r. The vector 
d denotes the position of a given molecule in the unit cell and is the number of the unit 
cells in the crystal. The Bragg scattering of a/j, in a large crystal is possible only when the 
momentum transfer is equal to one of the reciprocal lattice vectors 



k - k' = T . (23) 



The intensity of scattering peaks is determined by the Debye- Waller factor and by the value 
of (Jcoh for the specific momentum k and the scattering angle 'd. At high energies, the elastic 
Bragg scattering vanishes. Below the Bragg cutoff energy 

eB = lrlJMa^, (24) 

the condition (|23l) cannot possibly be fulfilled and the elastic coherent scattering disappears. 
The value of the shortest nonzero vector r is denoted here by Tmin- 

The Bragg-cutoff energies for various combinations of the muonic atoms and hydrogenic 
targets are shown in Table [Tll for the 3-K fee targets at zero-pressure. They have been 



calculated using the hydrogenic-crystal data from Refs. |22l, [23 
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TABLE II: The Bragg cutoff energy (in meV) for p/i, d(j,, and tjj, atom scattering in the 3-K fee 



crystals of nH2, nD2, and nT2 at zero pressure. 



lattice 


nH2 


nD2 


nT2 


Pli 


1.94 


2.13 


2.21 


dn 


1.02 


1.12 


1.17 




0.69 


0.76 


0.79 



When scattering takes place from a polycrystalline sample, which is usually the case in 
the muonic-hydrogen physics, Eq. (I^Tl) can be averaged over all orientations of the lattice. 
This gives the following cross section: 

X d \ L- —— - cosf 1 exp l-ZWiK^ 



5(^1 - ^ - cos^9^ exp [-2W{t 



The scattering now takes place in the Debye-Scherrer cones around the direction of k. The 
cones have the semi-angles subject to the condition cost? = 1 — r^/2A;^. In Fig. [3l ^9 is 




FIG. 3: (Color online) Scattering angle in Bragg scattering oi pfi in 3-K nH2 versus the collision 
energy and the value of the reciprocal-lattice vector r. The fee polycrystalline structure of the 
target is assumed. 



plotted as a function of the pfi kinetic energy in a 3-K nH2 target with the polycrystalline 
fee structure. The scattering angle is shown for the five smallest values r(z), z = 1, . . . , 5 of 
the vectors r. As one can see from Eq. and Fig. [31 the subsequent Bragg peaks appear 
at = 180°. When e increases, this backward scattering opens into a cone which moves 
continuously towards the forward direction. The total cross section cfohi^) expressed by 
the formula 

Att^ 1 ^^^'^ 1 
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where z{t) is the number of r- vectors with the same magnitude r and |-FVP denotes the 
average value of |F/v(t)P for just these vectors. 

In order to illustrate the energy dependence of interference effects, the function 



C 



Bragg 



\Fr\ 



, T<2k 

-Y- 



(27) 



is plotted in Fig. H] for p/i scattering in polycrystalline nH2 with the fee and hep struc- 
tures. This function is equal to the total Bragg cross section for scattering in a rigid lattice 
[exp(— 2iy) = 1], calculated per one molecule. It is assumed that cicoh is constant. The 
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FIG. 4: (Color online) The Bragg pattern in the case of scattering of a atom in a 3-K polycrys- 
talline nH2 target with the fee and hep structures. 



coherent scattering is forbidden below £:b ~ 2 meV. The subsequent Bragg peaks appear for 
the rising p/i energy. The Bragg-cutoff energy is slightly lower for the hep structure. 



V. COHERENT INELASTIC SCATTERING 

The inelastic coherent scattering is connected with the energy transfer between a/x and 
the coUective degrees of freedom of the target. Such a scattering describes the coherent 
annihilation and creation of phonons, but it does not include changes of the internal state 
of the target molecule (Ai? = 0). Thus, the energy transfer to the lattice is = e — e'. 
The inelastic coherent cross section for a/i scattering in a hydrogenic solid can be calculated 
using the methods developed for the coherent scattering of neutrons [2^. In this paper, 
coherent effects are taken into account only in the most important one-phonon processes. 
An estimation of the multiphonon cross sections is much more difficult. Therefore, these 
cross sections are calculated in the incoherent approximation. In the case of a Bravais 
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harmonic lattice, the coherent one-phonon cross section is given by the following formula: 



— E 



exp \-2W{k)\ J2 



C^coh 



X 



w 



(28) 



X n^^{q) 6(^uj + Wj) 6{k + q — t) 
+ (?^B(q) + l) <5('x^ - Wj) 6{k - q - r) 

q being the phonon momentum. There are three phonon-polarization vectors cr-'(q) {j = 
1,2,3) with the corresponding phonon energies Wj{q). The Bose factor for fixed u}j and q 
is denoted by n-'^{q)- The dispersion relation Wj = Wj{q) can be obtained by solving the 
lattice dynamics and is often measured in experiments, e.g., by means of neutron scattering. 
Because of the translational symmetry of the lattice, loj is a periodic function of q 

Wj{q)=Wjiq + T). (29) 

Thus, one can restrict values of q to the first Brillouin zone (r = 0) in order to know u}j at 
any q. The dispersion relation for small q takes the form 

Wj{q) = Cs{(T^) q , (30) 

where Cs{cr^) is the sound velocity in a given crystal. 

The cross section (I28p consists of the two terms. The first, which contains the expression 
S{uj + Wj)6{K + q — r) describes annihilation of one phonon. The term with the factor 
6{lj — Wj)5{n — q — t) corresponds to a/i scattering with simultaneous one-phonon creation. 
The annihilation processes vanish when the lattice temperature approaches zero, since there 
are no phonons at T = 0. The 5 functions in Eq. fl28l) represent conservation of both the 
energy and momentum, which is a basic feature of the coherent inelastic scattering. By 
virtue of Eqs. ([2]) and fl28l) . the initial and final energies and momenta of the muonic atom 
fulfill the conditions 

e' = e + Wj , k' = k + q — T , (31) 
in the case of one-phonon annihilation and obey the conditions 

e' = e — Wj , k' = k — q — r , (32) 

when one phonon is created. Therefore, for a fixed scattering angle, only phonons with 
a specific q and Wj{q) lead to the coherent one-phonon scattering. 

It can be shown that the coherent one-phonon creation near the forward direction k' ^ k 
(possible when r = 0) takes place only if the velocity of the impinging atom is greater than 
the sound velocity Cg [1^ . This means that the slowing down of the muonic atom via this 
process is impossible at the lowest energies. 

The isotropic Debye model of a solid is now used for the coherent phonon scattering. The 
average sound velocity Cg is substituted in Eq. fl30l) . which is taken as the first approximation 
to a real-crystal dispersion law. This assumption is reasonable because the coherent phonon 
processes are important in muonic atom deceleration only at the lowest energies. At higher 
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energies, incoherent scattering prevails. The factor \k ■ cr-?(q)p in Eq. (l28l) is replaced by its 
average value over a surface with a fixed w. Such an average for the cubic crystals equals 
which is a fair approximation even for most noncubic crystals [25|. Using the definition 

Z{w)^-^J2^[w-w,iq)] (33) 



of the normalized density of vibrational states, the summation over j and q in Eq. fl28l) is 
replaced by the integration 



3A^moi / dw Z{w) . 





The result of this integration is then averaged over the directions of r and q. Finally, one 
obtains the following formula: 



A^moi T c^coh exp(-2M/) 



X Z{UJ) ———Ujn{K,Uj) , 
iuj) + 1 



(34) 



where 



1 

TZ{k, u) = — — 0(tyD - |u;|) 
Vo Kq 

X ^ - 0(r - |k - q\) Q{k + q-T) 

T 

and G denotes the Heavyside function. There is a direct relation between this cross section 
and the incoherent one-phonon cross section 

(m£^L, = S:''<^-"K9^),./ '''' 

In the limit — 0, the function TZ is proportional to the factor 6{k, — r). This gives 
a geometrical condition similar to that for Bragg scattering in a polycrystalline target. Thus, 
in this limit, the coherent one-phonon scattering displays the same pattern of scattering 
peaks as that observed in the a/z-Bragg-scattering case (apart from r = 0). 
To describe the coherent phonon effects, the function 



V 



a 



inel 



is defined. A constant JF™"' is assumed here. In Figs. [5] and [6], both the coherent Ccoh and 
incoherent Cine functions are shown in the case of pfi scattering in the fee polycrystalline H2 
at T = 3 and 13.9 K. The one-phonon coherent and incoherent functions differ strongly 
below the Bragg-cutoff energy, especially at the lowest temperatures. On the other hand, 
the total coherent cross sections approach the corresponding incoherent ones at e greater 
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FIG. 5: (Color online) Functions Ccoh and Cine foi" scattering in the fee poly crystalline H2 at 3 K. 
Coherent one-phonon annihilation ("c— 1") and creation ("c+1") are shown using the solid lines. 
Incoherent (dashed lines) one-phonon processes have the labels: "i— 1" (annihilation) and "i+1" 
(creation). The labels "i— 2" and "i+2" stand for the two-phonon incoherent processes. 
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FIG. 6: (Color online) The same as in Fig. [Sj for T = 13.9 K. 



than a few meV. Since the energy spectrum of a created phonon is broad (^ 9 meV), the 
total cross section for one-phonon coherent creation is smooth. Oscillations of the coherent- 
annihilation cross section apparent at higher energies are due to a narrow width (^ 1 meV) of 
the energy spectrum of phonons which exist in the low-temperature targets. The functions 
Cinc(£), which describe the incoherent annihilation and creation of two phonons, are also 
plotted in Figs. [5] and [61 One can see that the incoherent approximation to multiphonon 
processes is reasonable because Ccoh for one-phonon coherent annihilation is much stronger 
than the correction from two-phonon incoherent annihilation. The only exception can be 
seen at e -C 1 meV in the 3-K target. However, this energy interval is not important when 
the slowing down of the p/x-atom from e ~ 1 eV is considered. 
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VI. EXAMPLES OF CROSS SECTIONS FOR HOMOGENEOUS HYDROGENIC 
TARGETS 



Some examples of the cross sections calculated for solid H2, D2, and T2 targets used 
in the TRIUMF, JINR, and RIKEN-RAL experiments are shown below. They include all 
the processes which preserve the total spin F of the muonic atom. In all the cases, the 
polycrystalline disordered (no specific orientations of molecules in the K = 1 state) fee 
structure has been assumed. At temperatures higher than 4 K, the targets can have the 
polycrystalline hep structure, which depends on the method of preparation and history of 



for the hep structure (apart from the Bragg cross section, which is calculated separately for 
the hep case). The reason for this is that the first three shells of neighbor molecules in 
both lattices are identical if the anisotropic interactions and orientations of molecules are 
neglected. The molar volumes of both structures are the same. The coherent phonon 
scattering is important below a few meV, where the acoustic-phonon approximation fl30l) is 
valid for the two structures, with practically the same value of the mean sound velocity. 

Although the Debye temperatures found in the literature are measured or calculated for 
crystals that do not exactly correspond to the targets used in muonic-hydrogen physics, 
these values are taken as a reasonable approximation. The inaccuracy involved by such an 
approach is estimated to be below 10%. When, for a specific condition, there is no data 
available, the Debye energy is calculated using the following equation: 



with the sound velocity Cg and the lattice constant Ofcc taken or calculated from Ref. [22|. 

In Figs. CHTSl the label "total" and the solid lines denote the total cross sections, which 
include all the coherent and incoherent processes. Contributions to the total cross sec- 
tion from the Bragg scattering (short-dashed line), phonon annihilation (dotted line, label 
"—phonon"), phonon creation (dash-dotted line, label +phonon"), and the rotational tran- 
sition K = 1 ^ (long-dashed line) are also plotted. The lines which describe the phonon 
processes do not include the rotational-vibrational transitions, although all the incoherent 
inelastic processes can take place with simultaneous phonon annihilation or creation. Thus, 
the rotational cross sections shown in these figures also contain the phonon terms. All the 
presented cross sections are calculated per one hydrogenic molecule bound in the solid tar- 
get. The cross sections for t/i atom scattering in zero-pressure solid nT2 at T = 3 and 20 K 
(just below the melting point) are shown in Figs. [7HTU1 For the ground state F = of the 
t/i spin, the scattering at the lowest energies is almost fully coherent (a small incoherence 
is caused by the slightly different single-molecule scattering amplitudes for the two lowest 
rotational states of T2). Below the Bragg cutoff, the elastic and phonon-creation coherent 
processes are forbidden. Thus, in this energy region, the coherent annihilation and incoher- 
ent phonon scattering are the only inelastic processes. At the lowest temperatures, however, 
these processes are strongly suppressed by the Bose factor. As a result, the total cross sec- 
tion in Fig. [7] falls by many orders of magnitude below ^b- This leads to a large increase of 
the an mean free path, which results in an enhanced emission of the cold (~ 1 meV) muonic 
atoms from the thin solid targets. This phenomenon has been already observed in the TRI- 
UMF experiment jll , in the case of p/i atoms. When the target temperature is raised, the 
low-energy phonon processes become more important and, near the melting temperature, 
the phonon cross section is quite large at e — (see Fig. [8]). 




However, the cross sections for the fee lattice are good approximations 




(37) 
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FIG. 7: (Color online) Total cross section for tfi{F = 0) scattering in zero-pressure solid nT2 at T = 
3 K. Contributions from Bragg scattering, phonon creation (-|-phonon), and phonon annihilation 
(— phonon) processes are also shown. 
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FIG. 8: (Color onhne) The same as in Fig. [7] for t^(F = 0) and T = 20 K. 



For F = 1, both the coherent and incoherent processes are significant since a-^^c is ap- 
preciable (cf . Table [T]) . Therefore, the total cross section below does not fall so greatly, 
compared to the F = case. The other important difference between the F = and F = 1 
cases is the rotational deexcitation K = 1 — >■ 0, which is due to the muon exchange between 
the tritium nuclei with the opposite spin projections. At the lowest energies, such a transi- 
tion leads to an effective acceleration of the muonic atom. This transition is forbidden when 
the total spin of tn{F = 0) is conserved in the collision process. On the other hand, such 
a transition can take place with the simultaneous excitation F = — 1 of the t/i spin state. 
However, the spin-flip process has a threshold of 0.237 eV (in the t^i + t center of mass) so 
that this reaction cannot occur at the lowest energies. 

Similar features of the cross sections appear in the pyU atom scattering in zero-pressure 
solid nH2. Figures [TT] and [T2] show the cross sections for F = and F = 1 at T = 13.9 K 
(below the melting ternperature of the solid H2). The analogous cross sections for T = 3 K 
were presented in Ref. js, 13]. Forp/i(F = 1), the incoherent scattering at the lowest energies 
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FIG. 9: (Color online) The same as in Fig.[7]for = 1) and T = 3 K. The label "1^0" denotes 
the rotational deexcitation K = 1 ^ oi the target molecule. 
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FIG. 10: (Color online) The same as in Fig. [9] for t/i(F = 1) and T = 20 K. 



is more pronounced than in the tritium case since Uinc is comparable to (Tcoh (cf. Table [T]). 

In the deuterium case, d/i atom scattering from the D2 molecule does not significantly 
depend on the ci/x-atom total spin 3]. Therefore, the cross sections for d[i{F = 1/2) and 
d[i{F = 3/2) scattering in solid D2 are very similar. The rotational deexcitation K = \ ^ ^ 
of the target molecule in d[i scattering which conserves the total spin F is allowed for both 
the F = 1/2 and F = 3/2 states. In Fig. [121 the total cross section for dii{F = 1/2) 
scattering in zero-pressure nD2 is plotted for the target temperature just below the melting 
point. For the ortho-D2 target, there is no rotational deexcitation. Thus, in this case, 
the phonon annihilation is the only (i/i-acceleration process, which is shown in Fig. [H] for 
T = 3 K. 

In Fig. [151 the total cross sections for dfi and scattering in 3-K solid nH2 are plotted. 
They can also be applied for solid H2 with a very small admixture of D2 or T2. When a fi~ 
beam is stopped in such a target, it is mainly pfi atoms which are created. Some of them can 
then hit a D2 or a T2 molecule, which mostly leads to the isotopic muon transfer [31|]. As 
a result, the released d^ or t/i atom has a kinetic energy on the order of a few tens eV. During 
the deceleration process, these atoms are emitted to vacuum while they reach the Ramsauer- 
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FIG. 11: (Color online) Total cross section for p^{F = 0) scattering in zero-pressure solid nH2 at 
T = 13.9 K. The notation is the same as in Fig. [7l 
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FIG. 12: (Color online) Total cross section for p^{F = 1) scattering in zero-pressure solid nll2 for 
F = at T = 13.9 K. The notation is the same as in Fig. [9l 



Townsend minima, which is apparent at about 10 eV. This mechanism has been used for 
the production of the energetic (1-10 eV) and t/i beams at TRIUMF [11, i, i, i, i, Q- 
The calculated differential cross sections establish the basis for accurate simulations of the 
emission of such muonic-atom beams. Fig. fTSl shows also that one should expect an enhanced 
emission of the cold (~ 1-meV) muonic atoms from the target, due to the falloff of the total 
cross sections below the Bragg cutoff. 



VII. MUONIC ATOM DECELERATION IN HYDROGENIC CRYSTALS 



When a /i~ beam is stopped in a hydrogenic target, the muonic hydrogen atoms are 
created in the highly excited Coulombic states. During the deexcitation process, these atoms 
gain kinetic energy. As a result, the ground-state IS atoms have a broad energy distribution 

10^ 



(see e.g., Refs. 



34, 35 



|36|), which can even extend to lO'* eV. Since the muonic 
atomic and muonic-molecular processes strongly depend on the a/i energy, it is necessary to 
know its time evolution. In the solid targets, the a^x deceleration from the highest energies 
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FIG. 13: (Color online) Total cross section for dfi{l/2) scattering in zero-pressure solid nD2 at 
T = 18.7 K. The notation is the same as in Fig. [9l 
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FIG. 14: (Color online) Total cross section for d/i(3/2) scattering in zero-pressure solid ortho-D2 
at T = 3 K. The notation is the same as in Fig. [71 



to about 0.1 eV is very fast (~ 1 ns). Therefore, only the last stage of afi slowing down, 
when solid-state effects are important, is considered here. The Maxwell distribution of the 
initial a/i kinetic energy with the mean value e^vg = 1 eV and the statistical population 
of the initial an spin states are assumed. The mean energy ^avg as a function of time and 
temperature has been evaluated using the calculated differential cross sections and Monte 
Carlo simulations. 

In Figs. [TBHT71 ^avg for the atom in solid nH2 and for the tfi atom in solid nT2 is plotted 
versus time. This energy is additionally averaged over the muonic-atom spin F. In fact, 
after about 1 ns, the higher spin states F = 1 of the pfj, and tfi atoms are depleted because of 
the large downwards spin- flip cross section [l6|. The deceleration of pfi^s to ^avg = 10 meV 



is very fast. This stage is much longer for t/z's as the cross section t/i(F = 0) + /!: is relatively 
small at the lowest energies [l3|- Below 10 meV, solid-state effects are very important and, 
therefore, they strongly affect the slowing down. This is especially visible at 3 K, when the 
phonon processes are suppressed and the classical thermal energy "ik^T /2 is much smaller 
than the Bragg cutoff energy (see Fig. [Tj). Below ^b, there is no effective deceleration 
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FIG. 15: (Color online) Total cross sections for d/x and i/x scattering in 3-K solid nH2. 
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FIG. 16: (Color online) Mean energy of a pfi atom in solid nH2 versus time at T = 3 and 13.9 K. 



mechanism. As a result, at the lowest temperatures, the mean energy of the muonic atoms 
in the steady state is much higher than the thermal energy. This energy is achieved after 
about 200 ns. The deceleration below 10 meV is slow since the scattering is dominated by 
the elastic Bragg process. On the other hand, when the target temperature tends to the 
melting point, 3k-BT/2 is greater than eb and the deceleration is faster. Also, the steady- 
state ^avg already approaches Sr (cf. the cross sections from Figs. [8] and [11]). One sees that 
the difference between the steady-state values of ^avg is much smaller than the difference 
of the corresponding classical thermal energies for the two limiting temperatures. This is 
especially pronounced in solid H2, where is the greatest and the melting temperature is 
the smallest, compared to the solid D2 and T2 targets. The analogous deceleration functions 
for atom in solid nD2 are plotted in Fig. [TSl Since the spin-flip rate for this case is very 
small [l6|, the deexcitation of the upper spin state is completed only after about 100 ns. This 
additionally extends the slowing down process as the hyperfine-splitting energy of 48.5 meV 
is released. 

In Fig. [191 the dependence of ^avg in the steady state is shown as a function of the target 
temperature, for the solid nD2 and nT2 targets. A full set of the differential cross sections 
has been calculated for about ten values of the temperature. Then, the mean energy has 
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FIG. 17: (Color online) Mean energy of a tfj, atom in solid nT2 versus time at T = 3 and 20 K. 
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FIG. 18: (Color online) Mean energy of a dfi atom in solid nD2 versus time at T = 3 and 18.7 K. 



been estimated by means of Monte-Carlo simulations. Below a few Kelvin, e^vg is almost 
constant, which is due to the lack of effective deceleration mechanisms at e < 1 meV. At 
higher temperatures, ^avg changes almost linearly and approaches 3/^6^/2 near the melting 
point. 



VIII. CONCLUSIONS 



The differential cross sections for muonic hydrogen atom scattering in the polycrystalline 
hydrogenic targets have been calculated using the scattering amplitudes for single molecules 
and the response function for the isotropic solid. A certain number of important approxi- 
mations has been made in the presented calculations. The scattering amplitudes for isolated 
hydrogenic molecules are obtained using the corresponding amplitudes for muonic atom 
scattering from free hydrogen-isotope nuclei [l^, 13] and the first-Born approximation. It is 
assumed that the molecular vibrations are harmonic and that there is no coupling between 
the vibrational and rotational degrees of freedom. The presence of electrons in the target 
molecules is taken into account by introduction of the effective electron-screening potential. 
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FIG. 19: (Color online) The calculated mean energy of a d/u atom in solid nD2 and of a t\i atom 
in solid nT2 at large times versus the target temperature. 



These approximations have been discussed in detail in Ref . 1^ . In the case of solid molecu- 
lar targets, it is assumed that the rotational-vibrational structure of the bound molecules is 



not changed, which is a fair approximation for the low-pressure hydrogenic crystals [2^, 123 
These targets are described in the harmonic approximation and the response function for 
the isotropic Debye solid is used. The Debye model is applied for estimation of the coher- 
ent and incoherent phonon cross sections. Coherent and incoherent effects are separately 
taken into account only in the most important one-phonon scattering. The cross sections 
for multiphonon processes are calculated in the incoherent approach. 

It happens that the muonic atom coherent scattering is very important in the H2, D2, 
and T2 targets at £ < 10 meV. In particular, a rapid decrease of the total cross section below 
the Bragg cutoff energy leads to an enhanced emission of cold muonic atoms at temperatures 
T < 5 K, where the phonon annihilation processes are strongly suppressed. The enhanced 
emission of p/z atoms has been observed and is well described by the calculated differential 
cross sections (sf. The coherent phonon creation disappears below about 2 meV. Since this 
process is the only mechanism of deceleration at the lowest energies, the mean energy of 
muonic atoms at T < 5 K is much greater than the corresponding energy in a perfect- 
gas target. Also, the deceleration process below 10 meV becomes much slower at such 
temperatures, since the elastic Bragg scattering in a heavy target cannot change the kinetic 
energy of the muonic atoms. 

Above about 10 meV, the rotational and then vibrational excitations of the target 
molecules support a very effective mechanism for the deceleration of muonic atom. At 
e > 100 meV, the cross sections for the hydrogenic crystals (per one molecule) tend to 
those calculated for the isolated hydrogenic molecules. However, the effective temperature 
of molecules in the hydrogenic crystals is much higher than the target temperature, due to 
the zero-point vibrations of the molecules in the lattice. This effect does not disappear (even 
at high collision energies). 

Although coherent and incoherent phonon effects have been estimated in the simplest 
approach of the isotropic Debye solid, the Monte Carlo simulations using the calculated 
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differential cross sections are in good agreement with the available experimental data. 
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